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A hybrid meson is a quark-antiquark pair in which, contrary to ordinary mesons, the gluon field 
is in an excited state. In the framework of constituent models, the interaction potential is assumed 
to be the energy of an excited string. An approximate, but accurate, analytical solution of the 
Schrodinger equation with such a potential is presented. When applied to hybrid charmonia and 
bottomonia, towers of states are predicted in which the masses are a linear function of a harmonic 
oscillator band number for the quark-antiquark pair. Such a formula could be a reliable guide for 
the experimental detection of heavy hybrid mesons. 
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I. INTRODUCTION 



The gluon being a colored object. Quantum Chromodynamics (QCD) allows the existence of exotic resonances, 
such as glueballs or hybrid mesons. Glueballs are composed of gluons only, while hybrid mesons contain a quark- 
antiquark pair as well. In particular, the study of heavy hybrid mesons is an active domain in theoretical and in 
experimental particle physics. BELLE and BABAR have already reported the discovery of several intriguing cc- or 
66-like resonances: One can quote the X(3872) [l],_but also the Y(4260) and T(10890) resonances iSl, that could 
be interpreted either as QQ hybrid mesons or as QQqq tetraquarks - see Refs. [3| for complete reviews. Notice that Q 
{q) denotes a heavy (light) quark. Moreover, one expects that future experiments like COMPASS, BESIII, GLUEX 
and PANDA should be very efficient in the detection of heavy hybrid mesons, especially of cc-type. 

There are two possible descriptions of hybrid mesons: First, a genuine three-body object made of a quark, an 
antiquark and a constituent gluon; Second, a two-body object made of a quark and an antiquark in the potential 
due to the gluon field in an excited state. In the framework of constituent models, it has been shown that these two 
pictures of the same object are, to a large extent, equivalent In this paper, we model the heavy hybrid meson 

as a QQ pair within an excited gluon field. 

In general, the string energy, and therefore the potential energy between the static quark and antiquark in the 
excited gluon field is given by H, ^ 

V{r) = ^Ja'^r'^ + 6, (1) 

where a is the usual string tension while 6 = iiraK + C is a term exhibiting the string excitation number K and 
a constant C . These values depend on the model adopted: a pure string theory Q or a more phenonienological 
approach 0,01 ■ In the present work, we choose the form predicted in Ref. 0, -f^ = 2ng -I- Ig and C = Stto, which is 
in very good agreement with lattice QCD for the standard value a = 0.2 GeV^ [13, [HI- Finally, 

6 = 27ra(2ng + /g + 3/2), (2) 

where Ug and Ig are respectively the radial and orbital quantum numbers of the constituent gluon simulating the 
excitation of the string. For the study of heavy hybrid mesons, it is therefore very interesting to calculate the 
eigenenergies of the Schrodinger equation governed by the potential ll]), or equivalently by the Hamiltonian 

H=^ + Va^r^ + b, (3) 
2/i 
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where = mQmQ/ {rriQ + mg) is the reduced mass and where the parameter h is given by Eq. ([2])- 

The aim of this report is to give an analytical expression for the mass of a heavy hybrid meson in the picture of 
an excited color field and to derive interesting physical consequences. Our analytical method relies on the auxiliary 
(or einbein) field method ( AFM) which has proved to be very powerful for such kinds of calculations [H, • The 
application of the AFM to Hamiltonian ([3]) is presented in Sec. |TT1 leading to an analytical mass formula. Using this 
formula, wc show in Sec. IIIII that it is possible to predict the general behavior of the heavy hybrid meson masses as a 
function of the quantum numbers of the system and search for possible towers of states. Our results are summarized 
in Sec. IIVI The excellent accuracy of the AFM mass formula is discussed in Appendix [XI through a comparison with 
exact numerical results. 



II. EIGENENERGIES 
A. Analytical expression 

Assuming that the potential ([1]) is the dominant interaction in a heavy hybrid meson, the mass of this system is 
■M'hin = mQ+rriQ+E, where E is an eigenvalue of the Hamiltonian ([3]). Using scaling laws (see Ref. [13]), dimensionless 
variables e and /3 can be defined as 

/o2\^/'^ 2/3 

e(/3) is an eigenvalue of the dimensionless Hamiltonian 



h=^ + y^M^. (5) 

Let us follow the general procedure of the AFM in order to find approximate expressions for the eigenvalues of 
the Hamiltonian ([5]). We first choose an auxiliary function P{x) — x^; the auxiliary field v is then defined by 

For the moment v is an operator, and Eq. ^ can be inverted to give cc as a function of x — I{v). Explicitly 



The AFM needs the definition of a Hamiltonian h{v) — + vP{x) + V{I{v)) — vP{I{v)). In our particular case, 

h{y) = ^ + iyx^ + — + i3v. (8) 

If we choose the auxiliary field in order to extremize h{i'), Sh/6iy\^=c — 0, then the value of this Hamiltonian for 
such an extremum is precisely the original Hamiltonian: h^v) = h. Instead of considering the auxiliary field as an 
operator, let us consider it as a real number. In this case, the eigenenergies of h are exactly known for all (n, I) 
quantum numbers: 

e{iy) = + ^ + (9) 

where, as usual, A^ = 2n + / + 3/2 is the principal quantum number of the state. 

The philosophy of the AFM is very similar to a mean field procedure. We first seek the value i^q of the auxiliary 
field which minimizes the energy, de{h') / dvliy^i^g , and consider that the value e(t'o) is a good approximation of the 
exact eigenvalue (T^ [isj . At this stage it is useful to define the new variable 

xo=N-'/%'^' (10) 

and the parameter 

Y = ^ (11) 

3^4/3- ^'') 
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The minimization condition is now concerned with xq and results from the fourth order reduced equation 

44 - 8a;o -3Y = 0. (12) 
The sohition of this equation can be obtained by standard algebraic techniques. It looks like 

xo - G(r) = ^^/V(Y) + l\/<ViY))-y^-ViY), (13) 

with 

/ ; \l/3 / , \-l/3 

V{Y) = [2+y^A + y3 j -Y [2+yU + Y^j . (14) 
Substituting this value into the expression of E{vo) leads to the analytical form of the searched eigenenergies, namely 

1 



G\Y) 



G{Y) 



(15) 



The problem is entirely solved. 

As it is shown in Ref. [l^, the same formula would be obtained for the choice P(r) = sgn(A)r'*' with A > —2, but 
with different forms for the quantity N . With the choice A = 2 made above, TV = 2n + / + 3/2. In this case, using 
results from Ref. [31 , it can be shown that formula (fT5)) gives an upper bound of the exact result. For A = — 1, 
N = n + l + 1 and the formula gives a lower bound. The qualities of these bounds are examined in Appendix |^ It is 
also shown in this appendix that the expression 

N ^ A{l3)n + l + C{l3), (16) 

with A{(3) and €{(3) defined by Eq. (|A2p . leads to an analytical formula which reproduces very accurately (up to 1%) 
the exact results. 



B. Asymptotic expansions 



Equation (jlSp is complicated but quite accurate. In order to get a better insight into this formula, it is interesting 
to compute several limits: 

^ N N'^ 

lim EAF = V/3 + —F^ — ; , (17) 

3iV2/3 p p2 

1.™/AF = + 22/3^,/3 - (18) 

The first two terms in the r.h.s. of Eq. P7|) are the solution of a harmonic potential, while Eq. with /3 = 
corresponds to the solution of a linear interaction obtained with the AFM [l3| . 

For the lowest excited state of the gluon field {K = 0) and for physical values of the parameters (a « 0.15-0.20 GeV^, 
rric « 1.1 GeV, ruh « 4.5 GeV), it comes that (3 « 7-20. In this case, (3 is large enough for the choice N = 2n-'rl-\- 3/2 
to be relevant (see Eq. (jA2p ). A harmonic oscillator band number B ^ 2n + I can thus be introduced to label the 
states. We can further assume that hybrid mesons with K — and B > A cannot be easily produced and discriminated 
from the low-lying K = I hybrids which would lie in the same mass range. So in the following, we will only consider 
that (3 £ [7, 20], B < 4, and K = 0. It is indeed more probable that hybrids with the lowest possible excitation of the 
gluon field will be first observed. 

As we are interested in the study of towers of states, it could be expected that the power expansion ([18]), valid for 
large N, should be used. But, /? is so large for heavy quarks that it is actually not the case. For example, when /3 — 20, 
one has Y ^ 5 ^ 1 even with the large value N = 10. By comparing with accurate numerical solutions obtained 



for the Hamiltonian h with the Lagrange mesh method [15| , it has been checked that the power expansion p?)) is far 



better. Moreover, the ratio of the third term over the second one in this expansion is at most 10%. So with very good 
approximation the eigenvalues of h are given by 

The approximate formula (llOp only depends on B. This exact degeneracy is actually broken by the non-harmonic 
character of the Hamiltonian h, but the breaking is small: It can be numerically checked that the maximal relative 
error of Eq. with respect to the exact values of e — ^/J3 is around 10%. So, B is a relevant classification number. 
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C. Final mass formula 



The Coulomb interaction must also play a role in heavy hybrids, since heavy quarks are expected to orbit close to 
each other and must feel strongly this short range interaction. The Hamiltonian ([S]) might then be completed with 
the potential nj x where 



, 2\ 1/3 
as lii"-^ 



6 V4a 



(20) 



as being the strong coupling constant and the 1/6 factor coming from color Casimir operator. Since as ^ 0.2 in the 
heavy meson sector 6], it appears that k < 0.02-0.07. For the large values of fi considered here, the contribution of 
the potential k/x can be computed in perturbation. Using the AFM results, one obtains 

^\ 2 « ^21) 



This contribution is at most around 4% of e — ^/]3 for cc hybrids and 10% for bh hybrids. So, in first approximation, 
it is relevant to neglect the Coulomb interaction in the search for towers of hybrid mesons. 
Under these conditions, the mass of a hybrid meson is given by 



Mhm ~ + mg + ^/2na{K + ^/T) 



''^^^ {B + 3/2). (22) 



MV27r(i^ + 3/2) 



For high excitation of the gluon field {K ^ 0), (3 increases and the accuracy of the approximation ([22]) improves. 
Using the AFM, an accurate mass formula for a ordinary meson with just the linear confinement is given by |12l | 

Q / 2 \ 1/3 / \ 



Afom~mQ + mQ + - (^-j \^—n + l + ^j . (23) 

Notice that = is not the ground state of the flux tube, but rather its first excited level. The ground state simply 
corresponds to 6 = in Eq. ([T|), that is a linear confining potential, and leads to formula 



III. PHENOMENOLOGY 



Some comments have to be done at this stage. Hamiltonian ^ describes a genuine heavy hybrid meson (no mixing 
with other hadronic states) within a spin-independent formalism. It will thus lead to qualitative global predictions 
rather than to a detailed mass spectrum. All spin effects are neglected but they should be weak for heavy hybrids 
since they are proportional to \/ {mQinQ) - see for example Ref. [l^ for a numerical check of that point. We think 
however that such global predictions are quite robust precisely because they do not depend on a fine-tuned model. 
Nevertheless, it is important to wonder whether our method will preferentially apply to some J^'-^ quantum numbers 
or not. Non-exotic quantum numbers, like 1 for example, must be examined very carefully because such quantum 
numbers allow for a possibly strong mixing with ordinary mesons. Such a mixing is by definition absent for exotic 
quantum numbers like 1 ^, 2^ , 3 ^ , . . . , although mixing with tetraquarks (or even glueballs) cannot be excluded. 
To our knowledge, the mixing between hybrid mesons and tetraquarks is far from being well-known theoretically, and 
we will not discuss it in the present work. 

Let us first estimate the value of the our model's parameters. The spin effects in mesons are minimal in the 1 
channel; it is thus relevant to fit the parameters by requiring Eq. (|23p to reproduce the radial trajectories of the 1 
cc and bh mesons. We first take a = 0.2 GeV^ so that potential (H]) optimally fits the corresponding lattice QCD 
data 0. Then, it is readily seen in Fig. [l]that the experimental data are well described by setting rac — 1.152 GeV 
and m.b — 4.620 GeV - the T(15') is poorly reproduced because of the neglect of the Coulomb term. Such values are 
not unrealistic when compared to the PDG values of 1.271q°^ GeV and 4.20j]QQy GeV [l^, and quite common in 
effective approaches. 

As already said in the introduction, the status of the 1^^ resonance F(4260) is not clear yet. If it is an ordinary 
meson, it should have S' = 1, n = 3 and I = [l^, and Eq. (|^ then leads to a mass of 4.512 GeV, quite far 
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Figure 1: Comparison between experimental cc and hh radial trajectories (full circles) and the predictions of Eq. (|23|l for 
a = 0.2 GeV^, rric = 1.152 GeV, and vrih = 4.620 GeV. Experimental data come from the PDG [l^ and the quantum numbers 
assignment of Ref. [l3| for the cc mesons is followed. Dashed lines are plotted to guide the eyes. 



from the experimental value. If it is a hybrid meson, the quantum numbers must be S' = 1, n = and I = 1 [20| . 
and Eq. ([22|) leads to a more similar mass of 4.239 GeV. But, the repulsive Coulomb term would slightly decrease 
this agreement. The ordinary meson interpretation for the 1 resonance T(10890) has been discarded in Ref. [1]; 
Eq. ([22]) is compatible with the hybrid meson picture since it leads to a mass of 10.894 GeV for the corresponding 
state. Our model does not disagree with a hybrid meson assignment for both states, but no definitive conclusion can 
be drawn at this stage, mainly because of the neglect of the mixing with ordinary mesons. 

We now propose a modus operandi to separate the hybrid mesons from the ordinary ones. Let us first assume that 
we have a sufficient amount of experimentally found QQ-like resonances - this should be possible in a near future 
thanks to the forthcoming experiments we mentioned in the introduction. Then, for each J^*-^ state, one has to find 
the lowest possible value of B that is compatible with these quantum numbers in a hybrid meson picture. We recall 
that P = (— 1)'+'9 and C = (—iy+^iq+^ in this case, and that Ig — when K = 0. Once this step is achieved, one 
can plot the experimental masses versus B. Our prediction is that the hybrid mesons, or at least states that are 
dominated by a hybrid meson component, will be located along a straight line while other states like ordinary mesons 
or tetraquarks will not (typically for example, one has Mom — — "m-Q oc n^^^ and P^^). 

Formula (l22ll can be rewritten as follows for A' = 



Mhm(B) = W (B - Bo) + Mhm(So). (24) 

y hvStt 

It is plausible that an exotic 1 ^ hybrid meson, thus with B = 1, will be first discovered. Effective approaches as 
well as lattice QCD indeed reach the conclusion that the 1 ^ state is the lightest exotic QQ hybrid meson (see e.g. 
Refs. [H, [2l|). Then, by setting Bq = 1 and Mhm(So) = Mi-+, Eq. pi)) states that the first exotic states will be 
located on a straight line, as schematically represented in Fig. [21 Note that in our ap pro ach, the 1 ^ cc state lies 
below the experimental estimation of the DD* threshold, which is of about 4.3 GeV {D* denotes the P-wave 
excitation of the D meson). 

The results of formula (|24p can also be compared with lattice QCD. Most of the efforts in lattice QCD were 
devoted to the computation of the lowest-lying 1 ^ hybrid meson masses. We can quote: 4.420±0.013 GeV [U 



4.369 ± 0.136 GeV 123|, and 4.405 ± 0.038 GeV |24| for the 1"+ cc states, and 10.82 ± 0.08 GeV [25| and 10.977 ± 
0.123 GeV [l^ for the 1 ^ bb states. It is worth saying that the 1 and 1~~ bb hybrid mesons are degenerate in 
Ref. [1^, showing the weakness of spin effects for systems of bottom quarks. Formula is in agreement with the 
current estimations as shown in Fig.[2l it would be very interesting that other masses become available in lattice QCD 
in order to check whether they are located on the straight lines we predict or not. 
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Figure 2: Masses of exotic QQ hybrids versus the band number B as predicted by formula (|24|l (circles and dashed lines). 
Parameters of Fig.[T]are used. An explicit example of exotic J^*^ trajectory is given; the lattice data of Refs. [2^. l23l. [23. [25l. l26l] 
are also plotted for comparison (crosses). 



IV. CONCLUSION 



Starting from a potential model of hybrid mesons which is in agreement with lattice QCD, we predict that the mass 
of the lowest-lying genuine heavy hybrid mesons is a linear function of the band number B — 2n + I, n and I being 
the quantum numbers of the quark-antiquark pair. The slope and the intercept at the origin depend on the excitation 
quantum number K of the gluon field. These heavy hybrid mesons form towers of states organized as lines in plots 
where the masses are presented versus the harmonic oscillator band number. We believe that this property can be 
an interesting tool to disentangle resonances which are dominated by a heavy hybrid meson component from other 
hadronic states in future experiments, just as the existence of Regge trajectories is an important guide to identify light 
mesons. Apart from experiment, an important check of the present results would be the computation of higher-spin 
exotic hybrid mesons on the lattice. We hope that such results will be available in the future. 
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Appendix A: COMPARISON TO EXACT RESULTS 



The AFM cannot give strong constraints on the dependence of N in terms of {n,l). In particular, had we chosen 
P(r) = sgn(A) r^, the better choice for N would have been TV = A{X)n + I + C(A), with the quantities A{\) and C(A) 
given in Ref. (l^ . The square root potential y/ x'^ + (3 ensures a smooth transition from a linear form (A = 1 but in 
this case we have only approximate expressions) to a quadratic form (A = 2 and in this case the values are exact) as 
(3 increases from to oo. 

The procedure we adopt is based on the following points: 

• We calculate the exact values e*;(/3) for < n < nmax, < Z < l^^x and for a given set of (3 values. This 
program is achieved using a very powerful method known as the Lagrange mesh method (described in detail in 
Ref. [IBl)- For our purpose, we consider that rtmax — ^max = 4 is a good choice. For any calculated value, we 
have an accuracy better than 10~^. 

• We calculate the approximate values £„;(/?) using Eqs. with Y given by Eq. pT|) in which TV is deduced 
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from Eq. ITBl) . for the same set of j3 values. Building the x-square 



we request the coefficients A{f3) and C(/3) of iV to minimize this function. The obtained values are represented 
by black dots in Figs. [31 

• In order to obtain functions which are as simple as possible, continuous in /5, and which reproduce at best the 
above calculated values, we choose hyperbolic forms and require a best fit on the set of the sample. Explicitly, 
we find 

8^+ 102 30/^ + 53 

4/3 + 57' 20/3 + 39- 

These integers are rounded numbers whose magnitude is chosen in order to not exceed too much 100. The 
corresponding values are plotted as continuous curves in Figs. [31 They have been constrained to exhibit the 
right behavior A ^ 2 and C 3/2 for very large values of (3. Formulas ((X2|) give A(0) = 102/57 w 1.789 
and C(0) = 53/39 « 1.359. These values are such that A{0) « tt/VS w 1.814 and C(0) « y37r/4 « 1.360, as 
expected from the results of Ref. [12] in the case of a nonrelativistic linear potential. 



Table I: Comparison between the exact values eJ^i(/3) (2nd line) and analytical approximate expressions e-niiP) for the eigenvalues 
of Hamiltonian ^ with f3 — 1. For each set (n, I), the exact result is obtained by numerical integration. 3rd line: approximate 
results are given by Eqs. (|15|l with Eqs. (Illf) . p6p and (IA2P : 1st line: upper bounds obtained with N = 2n + I + 3/2; 4th line: 
lower bounds obtained with N — n + I + 1. 





I = 


I = 1 


I = 2 


; = 3 


I = 4 


n = 


1.94926 


2.49495 


2.99541 


3.46197 


3.90193 




1.91247 


2.45074 


2.94841 


3.41419 


3.85430 




1.89549 


2.44621 


2.95032 


3.41969 


3.86189 




1.65395 


2.22870 


2.75000 


3.23240 


3.68492 


n = 1 


2.99541 


3.46197 


3.90193 


4.32027 


4.72059 




2.89556 


3.34652 


3.77899 


4.19405 


4.59335 




2.85420 


3.32970 


3.77678 


4.20097 


4.60620 




2.22870 


2.75000 


3.23240 


3.68492 


4.11355 


n = 2 


3.90193 


4.32027 


4.72059 


5.10556 


5.47723 




3.74112 


4.14232 


4.53310 


4.91307 


5.28251 




3.69078 


4.11913 


4.52783 


4.91998 


5.29790 




2.75000 


3.23240 


3.68492 


4.11355 


4.52250 


n — 3 


4.72059 


5.10556 


5.47723 


5.83725 


6.18692 




4.50374 


4.87138 


5.23246 


5.58628 


5.93264 




4.44883 


4.84403 


5.22459 


5.59242 


5.94903 




3.23240 


3.68492 


4.11355 


4.52250 


4.91485 


n = 4 


5.47723 


5.83725 


6.18692 


6.52732 


6.85935 




5.20859 


5.55148 


5.88996 


6.22329 


6.55111 




5.15078 


5.52098 


5.87970 


6.22821 


6.56756 




3.68492 


4.11355 


4.52250 


4.91485 


5.29295 



Since our results are exact for f3 —^ oo, one has obviously x = in this limit. The error is maximal for small 
values of /? but, over the whole range of /? values, the results given by our analytical expression can be considered as 
excellent. Just to exhibit a quantitative comparison, we report in Table [J the exact e*;(/3) and approximate eni{/3) 
values obtained for f3 — 1, a, value for which the corresponding potential is neither well approximated by a linear 
one nor a harmonic one. As can be seen, our approximate expressions are better than 1% for any value of n and I 
quantum numbers. Such a good description is general and valid whatever the parameter /3 chosen. 

The upper bounds obtained with P(r) = are far better than the lower bounds computed with P{r) = — l/r. This 
is expected since the potential \/ a'^r^ + b is closer to a harmonic interaction than to a Coulomb one. Better lower 
bounds could be obtained with P{r) = r. But, the exact form of N is not known for this potential, except for I = for 
which TV can be expressed in term of zeros of the Airy function. With the approximate form N = {n / \/3)n + 1 + ^/Sn / 4 
[H, 111] , we have checked that results obtained are good but the variational character cannot be guaranteed. 
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Figure 3: Best values of the coefficients A{I3) and C(/3) to parameterize the eigenvalues of Hamiltonian ([5]): numerical fit with 
Eq. (|XT|) (dots); functions (E2)) (solid line). 
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